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Abstract. We prove Lieb-Robinson bounds and the existence of the thermo- 
dynamic limit for a general class of irreversible dynamics for quantum lattice 
systems with time-dependent generators that satisfy a suitable decay condition 
in space. 

1. Introduction 

For a quantum many-body Hamiltonian describing bulk matter we expect that 
the Heisenberg dynamics converges in the thermodynamic limit to a well-defined 
one-parameter flow of transformations on the observable algebra. Early results of 
this kind were obtained for quantum spin systems j31L 1261 128] , which were followed 
by generalizations that included examples of irreversible dynamics described by a 
semigroup of completely positive unit preserving maps [9l 130] . See, e.g., [61 [8l 
[3l 1321 I17j for results on the thermodynamic limit of a number of examples of 
semigroups of completely positive maps. In this work we study a general class of 
irreversible dynamics for quantum lattice systems with generators that are sums of 
bounded terms that may depend on time and that satisfy a suitable decay condition 
in space. 

Following the argument of [27] propagation bounds of Lieb-Robinson type |15| 
have recently been used to prove a number of new results on the existence of 
the thermodynamic limit [181 \T\ I20L [2- These recent developments were made 
possible by extensions and improvements of the Lieb-Robinson bounds themselves 
[TTl[2n[T3l[T9l [22]. 

Lieb-Robinson type bounds for irreversible dynamics were, to our knowledge, 
first considered in |12j in the classical context and in |25j for a class of quantum 
lattice systems with finite-range interactions. Here, we will extend those results by 
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proving a Lieb-Robinson bound for lattice models with a dynamics generated by 
both Hamiltonian and dissipative interactions with suitably fast decay in space and 
that may depend on time. See Assumption [1] for the precise conditions. Then, we 
use our result to prove the existence of the thermodynamic limit of the dynamics in 
the sense of a strongly continuous one-parameter flow of completely positive unit 
preserving maps. 

Our results are applicable to a wide range of model systems in statistical me- 
chanics, quantum optics, and quantum information and computation. In each of 
those areas, it is often necessary to incorporate dissipative and time-dependent 
terms in the generator of the dynamics. Fortunately, there is a large number of 
interesting systems defined on a lattice, which so far is the only setting accessible 
by our methods to prove Lieb-Robinson bounds. It is probably not a coincidence 
that proofs of the existence of the thermodynamic limit of the dynamics have so far 
also been mostly restricted to lattice systems. Here, 'lattice' has to be interpreted 
loosely to mean a discrete set of points that are typically thought of as distributed 
in space. In the case of the positions of atoms in a crystal, these positions can in- 
deed be described by a lattice, but all one needs is the structure of a metric graph 
satisfying some regularity conditions. The detailed setup is given in Section El 

The existence of the thermodynamic limit is important as a fundamental prop- 
erty of any model meant to describe properties of bulk matter. In particular, 
such properties should be essentially independent of the size of the system which, 
of course, in any experimental setup will be finite. In the past five years, Lieb- 
Robinson bounds have been used to prove a variety of interesting results about 
condensed matter systems. See j23j for a brief overview of the applications of 
Lieb-Robinson bounds. 

The paper is organized as follows. First, we describe the general setup necessary 
to state our main results, which we do in Section [2l In that section we also state 
the three main theorems we prove in this paper. Theorem [T] states that solution of 
the differential equation (master equation) defined by finite volume generators we 
consider is a well-defined quantum dynamics, i.e., a continuous family of completely 
positive unit preserving maps on the algebra of observables. The proof of this 
theorem is obtained by standard methods, but for completeness we included it here 
in Section[3l Theorem[2]is the Lieb-Robinson bound, i.e., the propagation estimate 
for irreversible dynamics. Again the theorem is stated in Section[2]and then proved 
in Sectional Theorem |31 the existence of the thermodynamic limit, is proved in 
Section [51 



2. Setup and main results 

We consider quantum systems consisting of components associated with the 
vertices x € T, where F is a countable set equipped with a metric d. We assume 
that there exists a non-increasing function F : [0, oo) — > (0, oo) such that: 
i) F is uniformly integrable over F, i.e.. 



:= sup Y.F{d{x,y)) < oo, 
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and 

ii) F satisfies 

^ Fidix,z))F{d{y,z)) ^ 
C := sup } — — 7- < oo. 

Having such a set F and a function F that satisfies i) and ii) , we can define for 
any /j, > the function 

F^{x) = e-^-F{x), 

which then also satisfies i) and ii) with ||F^|| < ||F|| and < C. 

The Hilbert space of states of the subsystem at a; S F is denoted by T-Lx- For 
any finite subset A C F the Hilbert associated with A is 

The algebra of observables supported in A is defined by 

AA = (^B{nx), 

where B{Hx) is the set of bounded linear operators on 'Hx- If Ai C A2, then we 
may identify Aai in a natural way with the subalgebra Aai ® ]1a2\Ai of Aa2^ and 
simply write Aai C Aa2 ■ The algebra of local observables is then defined as 

-^r' =[_]Aa. 

Acr 

The C*-algebra of quasi- local observables is the norm completion of ^j?'^. See 
[4t [5] for more details about this mathematical framework. 

The support of the observable A G is the minimal set AT C A for which 
A = A' ® tA\x for some A' & Ax- 

The generator of the dynamics is defined for each finite volume A C F, and, 
in general, contains both Hamiltonian interactions and dissipative terms, which we 
allow to be time-dependent. The Hamiltonian terms are described by an interaction 
i>(f, •) which, for all t S R, is a map from a set of subsets of F to At, such that 
for each finite set X CT, <^{t,X) G Ax and $(t,X)* = <^{t,X). The dissipative 
part is described by terms of Lindblad form determined, for each finite A" C F, 
by a set of operators La{t, X) G Ax-, a = 1, • • ■ , ^{X). We can allow for the case 
N(X) = 00, if we impose a suitable convergence condition on the resulting series 
for the generator. Then, for any finite set A C F and time i G M we define the 
family of bounded linear maps Ca ■ Aa Aa, i-C. Ca G B{Aa, Aa), as follows: 
for aU A e Aa, 

(2.1) ^zmA) = tm,z),A] 



N(Z) 

J2 [Ll{t,Z)ALa{t,Z) - -{La{t,zyLa{t,Z),A]] 



(2.2) CA{t){A) = ^z{i){A), 



ZdA 



where {A, _B} — AB + BA, is the anticommutator of A and B. The operators 
^z{t) can be regarded as bounded linear transformations on Ax, for any AT C A 
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that contains Z, which arc then of the form 'i>z{t) ® id^^^^- The norm of these 
maps, in general, depends on X , but they are uniformly bounded as follows: 

NiZ) 

||*^(t)||<2||$(t,Z)|| + 2^ \\Ut,Z)r 

a=l 

If N{z) = oo, we can insure uniform boundedness by assuming that the sums 
converge. It is more general and more natural, however, to as- 
sume that the maps ^'^(i), defined on Az are completely bounded. By definition, 
^' G B{Az) is called completely bounded if for all n > 1, the linear maps 4' ® idM„, 
defined on Az ® M„, where Af„ ~ i3(C") are the n x n complex matrices, are 
bounded with uniformly bounded norm. This means that we can define the cb- 
norm of 5* by 

||1'||cb = sup||*(»idMj| <oo 
ji>i 

In particular, this definition implies that the cb-norm of ^z{t), which can be con- 
sidered as a linear map defined on for all A C F such that Z C A, is independent 
of A. See [7 tl24j for more information on completely bounded maps. Assuming that 
||^z(i)||cb is finite is more general than assuming that the series J2T=i \\^a{t, Z)\\^ 
converges which, however, is a useful sufficient condition for it. In particular, there 
are situations where the sum in (|2.ip only converges in the strong operator topology 
but nevertheless yields a well-defined limit with finite cb-norm. 

Assumption 1. Given {T,d) and F as described at the beginning of this section, 
the following hypotheses hold: 

(1) For all finite A C F, C\{t) is norm- continuous in t, and hence uniformly 
continuous on compact intervals. 

(2) There exists fi > such that for every t G R 

(2.3) 11*11,, := sup sup J2 < ^■ 

where \\ ■ ||cb denotes the cb-norm of completely bounded maps [24j . 
Note that 

l|/:A(i)ll < E ^ E E ll*zWllcb<ll*llt,,|A|||F|j. 

Z<ZA x,yeAZ3x,y 

We define 

(2.4) M, = ||f||,,,|A|||J^||. 

Then by ([^ one gets < Mt for s <t. 

Fix T > and, for all A e Aa, let A{t),t G [0,r] be a solution of the initial 
value problem 

(2.5) f^A{t) = CA{t)A{t), AiO)=A. 

Since ||£A(i)|| < Mt < oo, this solution exists and is unique by the standard 
existence and uniqueness results for ordinary differential equations. For < s < 
t <T, define the family of maps {7t^s}o<s<t C B{Aa, Aa) by 7^^(A) = A{t), where 
A{t) is the unique solution of ()2.5|) for t G [s,T] with initial condition A{s) = A. 
Then, the cocycle property, ^t,s{A{sj) ~ A{t), follows from the uniqueness of the 
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solution of p.Sp . Recall that a linear map 7 : ^ — > S, where A and B are C*- 
algebras is called completely positive if the maps 7 (g) id : ^ (81 M„ — > S M„ 
are positive for all n > 1. Here M„ stands for the n x n matrices with complex 
entries, and positive means that positive elements (i.e., elements of the form A* A) 
are mapped into positive elements. See, e.g., |24| for a discussion of the basic 
properties of completely positive maps. In particular, we shall use the property 
that every unit preserving (i.e. 7(1^) = Ig) completely positive map 7, is a 
contraction: ||7(A)|| < ||A||. 

As a preliminary result we prove the following Theorem [1] in Section [31 It 
extends the well-known result for time-independent generators of Lindblad form 
[16j to the time-dependent case. 

Theorem 1. Let A be a C* -algebra, T > 0, and for t e [0,T], let C{t) be a norm- 
continuous family of bounded linear operators on A. If 

(i) C{t){t)^Q- 

(ii) for all A e A, C{t){A*) = C{t){A)* ; 

(Hi) for all AeA, Cit){A*A) - C{t){A*)A ~ A*C{t){A) > 0; 

then the maps jt,s, < s < < < T, defined by equation 112. 5\) . are a norm- continuous 
cocycle of unit preserving completely positive maps. 

It is straightforward to check that the C\{t) defined in ()2.2p satisfy properties 
(i) and (ii) . Property (iii) , which is called complete dissipativity, follows immediately 
from the observation 

N{Z) 

C{t)iA*A)-C{t){A*)A-A*C{t){A)=J2 E t^' ^-(*' ^)nA ia(i, ^)] >0. 

ZCA a=l 

Therefore, using this result, we conclude that, under Assumption [U for all finite 
A C r, the maps 7^^, < s < t, form a norm-continuous cocycle of completely 
positive and unit preserving maps. 

Section U is devoted to proving a Lieb-Robinson bound for the irreversible 
dynamics 7^^. For reversible dynamics given by the one-parameter group of au- 
tomorphisms Tt describing the Heisenberg dynamics generated by a Hamiltonian, 
Lieb-Robinson bounds take the following form: there are constants v, n > such 
that for A e Ax and B £ Ay, 

(2.6) \\[A,Tt{B)]\\ < C(A,B)e-'^('^(^'^)-"'*l) , 

where d{X,Y) denotes the distance between X and Y and C{A,B) is a prefactor, 
which typically has the form c\\A\\ \\B\\ min(|Ar|, |F|), for a suitable norm [| • [| on 
the observables A and B, and a suitable measure | • | on the size of the supports X 
and Y. Bounds of this form are sufficient to determine the approximate support of 
the time-evolved observable Tt{B). See, e.g., [2j Lemma 3.1]. 

For irreversible dynamics, it turns out to be both natural and convenient to 
consider a slightly more general formulation. For X C A, let Bx denote the sub- 
space of B{Ax) consisting of all completely bounded linear maps that vanish on 1. 
See the discussion directly preceding Assumption [T] for the definition of complete 
boundedness and the cb-norm || • |[cb. It is important for us that all operators of 
the form 

KxiB) - [A, B] +Y,[L:BLa - -{KLa, B}] , 
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where A, La G Ax, belong to Bx, with 

N 

\\JCxU<2\\A\\+2Y,\\Laf. 

a=l 

In particular, operators of the form [A, ■] appearing in the standard Lieb-Robinson 
bound (|2.6[) are a special case of this general form. Then, we can regard ICx as 
a linear transformation on Az, for all Z such that X G Z, hy tensoring it with 
id^^^^-, and all these maps will be bounded with norm less then H/Cx ||cb- 

Theorem 2. Suppose Assumption[l\ holds. Then the maps •y^^ satisfy the following 
bound. For X,Y d A, and any operators fC G Bx o,nd B g Ay we have that 

Note that the bound in this theorem is uniform in A. This is important for the 
proof of existence of the thermodynamic limit of the dynamics, which is the main 
application of Lieb-robinson bounds in the present paper. 

As a final comment about the use of the cb-norm in the definition of ||'I'||t,p (see 
(|2.3p ) we would like to point out that volume-independent bounds for the operator 
norm, such as ||[^', -JH < 2||\1/||, which appear in all previous Lieb-Robinson bounds, 
are always an upper bound for the norm used here. This is also true for the case of 
reversible dynamics. The bound stated here will give a sharper result in some cases. 
It has been suggested that the addition of dissipative terms to the generator of the 
dynamics would not increase the Lieb-Robinson velocity. For example, in j29| it is 
correctly argued that the bounds derived in that paper remain valid without change 
if one adds particle loss terms. While this is trivially true for arbitrary single-site 
terms, it is not clear that the same comparison would hold in general. We also need 
to alert the reader that the bounds derived in [29], while valid for lattice bosons 
with a finite number of particles, have a prefactor which depends on the particle 
number (at least linearly for the simplest observablcs, and worse than linear for 
more general observablcs). In this sense the results of |29j are not a true extension 
of [25j to many-body boson systems as discussed, e.g., in the recent book [33j . 

Also note that the bound given in Theorem [2] can be further improved by 
omitting in the definition of ||\E'||t,^ all terms that act on a single site, i.e., which 
belong to B^^j for some x G F, and also all terms in Bx, where X is the set for which 
the bound is derived. By the argument in |19| we can even allow the single site 
terms to be unbounded, as long as they lead to a well-defined single-site dynamics. 

In this paper, we restrict ourselves to applying the Lieb-Robinson bound of 
Theorem [5] to proving the existence of the thermodynamic limit of a general family 
of irreversible dynamics. Other applications, such as approximate factorization of 
invariant states, analogous to what is done for ground states of reversible systems 
in [To] . 

The setup for the analysis of the thermodynamic limit can be formulated as 
follows. Let F be an infinite set such as, e.g., the hypercubic lattice Z'^. We prove 
the existence of the thermodynamic limit for an increasing and exhausting sequence 
of finite subsets A„ C F, n > 1, by showing that for each A G Ax, (7t^r(^))">i 
is a Cauchy sequence in the norm of ^r- To this end we have to suppose that 
Assumption [1] (2) holds uniformly for all A„, i.e., we can replace A in p.3|) by F. 
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Theorem 3. Suppose that Assumption]^ holds and, in addition, that i2. 3\) holds 
for A = r. Then, there exists a strongly continuous cocycle of unit-preserving 
completely positive maps ^ on Ar such that for all < s < t, and any increasing 
exhausting sequence of finite subsets A„ C T, we have 

(2.7) lim ||7t"(^)-7M(^)ll=0> 

n— >oo 

for all A e Ar- 

3. Finite volume dynamics 

Let C{t), t > 0, denote a family of operators on a C*-algebra A satisfying the 
assumptions of Theorem [T] and for < s < i consider the maps A 3 A i-^ 'yt,s{A) 
defined by the solutions of (|2.5I) with initial condition A at t = s. Without loss of 
generality we can assume s = in the proof of the theorem because, if we denote 
C{t) = C{t + s), then jt.s ~ 7t-s,o, where 74,0 is the maps determined by the 
generators 

The maps 3 satisfy the equation 

(3.1) 7t,s=id+/ C{T)'-fr,sdT. 

J s 

In our proof of the complete positivity of 74.0 we will use an expression for 74^0 a-s 
the limit of an Euler product, i.e.. approximations T„{t) defined by 

(3.2) T„(t) = n fid+i/:(^: 

k=n ^ 

The product is taken in the order so that the factor with = 1 is on the right. 

Lemma 1. Let C{t), t > 0, denote a family of operators on a C* -algebra A satis- 
fying the assumptions of Theorem]^ Then, uniformly for all t G [0,T], 

lim ||T„(0-74,o|| -0, 



where T„(t) is defined by I13.S\) . 

Proof. From the cocycle property established in Section [51 we have 

1 

7t,Q = n ^t#,t^- 

k—n 

Now, consider the difference 

1 1 

Tn{t) - 7*.0 = n (id+^/:(^)) - n 74|,4i^ 
k—n k—n 

n rj-\-i 



E 



n(id+^/:(t^)) 



.k= 



id+^/:(i^))-74, 



7tiiii,o- 



To estimate the norm of this difference we look at each factor separately. 

Using the boundedness of C{t) and the fact that Mt, defined in ()2.4|) . is in- 
creasing in t, the norm of the first factor is bounded from above by 

II 'fiiid+Lc{t^m < n(i + -\\cit^)\\) < (1 + '-M,r. 

k—n k—n 
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To bound the second factor notice that from (13.11) we obtain 



||7t,.||<l + ||/:(r)||||7t,s|Mr. 

Then by Gronwall inequahty |14[ Theorem 2.25] we have the following bound for 
the norm of the jt,s' 

Using again p.ip we can rewrite the second factor as follows: 



n n 



n n ' 



(7ti,ti-i - id) 



3-1 



ds 



C{f^—) - C{s) Us - / C{s) / £(T)7^,3^drds 



Therefore, the second factor is bounded from above by 



n n 



t , 



3-1 /iJ-1 



3-1 



2n / ' 



where e„ ^-Oast/n— >0 due to the uniform continuity of C{t) on the interval [0, t]. 
The third factor can be estimated in a similar way: 



ii7t2^,oii = n iiTti.t 



k=j-i 
1 



^ ^ n n 



< n ii + -\\c{s,i-m 



< (1 + -Mt)". 
Therefore, combining all these estimates we obtain 

\\Tnit) - jt,o\\ < ^(1 + -MtT- fe„ + M!e^''^^] (1 + -MtT 
n n \ 2n I n 



2n 



This bound vanishes as n — > c». 



□ 
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To prove Tlieorem[l]we use the Euler-type approximation established in Lemma 
[H We show that the action of T„ on a positive operator gives a sequence of bounded 
from below operators such that the negative bounds vanish as n goes to oo. 

Proof of Theorem [T] First, we look at the each term in the Euler approx- 
imation r„ separately. For any t and s the complete dissipativity property (iii) of 
C{s), assumed in the statement of the theorem, implies 

< {id+tCis)){A*){id+tCis)){A) = {A* +tCis){A*)){A + tCis){A)) 
= A* A + tA*C{s){A) + tCis){A*)A + fC{s){A*)Cis){A) 
< A* A + tCis){A*A) + t^C{s){A*)Cis){A). 

Since {C{s){A)y {C{s){A)) < \\Cis)f\\A\\, one gets 

(3.3) < iid+tCis)){A*A) + e-\\C{s)f\\Af 

(3.4) < {id+tCis))iA*A) + fMf\\Af. 

Let us apply the above inequality to the operator^, where B*B := j|j4|p — 
Note that \\B*B\\ < \\A\\^, so ||B|| < 

(3.5) < {id +tC{s)){\\Af - A* A) +t^M^\\Af 

(3.6) = II Af - (id +tC{s)){A* A) + fM^WA^ 

From the p.3p and p.Sp we obtain 

(3.7) -fM^\\A\\^ < {id+tCis)){A*A) < {l + t^M^)\\Af 
and therefore: 

-{l + t^M^)\\Af < {id +tC{s)){A* A) < {l + t^M^)\\Af. 

So we get 

(3.8) \\iid+tC{s)){A*A)\\<{l + tHl^)\\Ar . 

Now, in order to bound the approximation T„ we first derive the following 
auxiliary estimate. For any fixed n > I we have: 

1 ^ n-l 

(3.9) l[{id+sC{ks)){A*A) > ^s^ArMl{l + 

k—n k—0 

where the value of 5 is chosen to be such that 

(3.10) D{s) := 1 + s'Ml < (1 + ^)"- 7(1 + ^)""', 

n — 1 n ~ 2 

with the convention that (1 H — ^)"^^ — 1, for n = 1. 

We prove this claim by induction. The statement holds for n = 1 by p.Sp . 
Now, assume that (|3.9p holds for n — 1. Then 



1 ^ n-2 

n (id ^sL{ks)){A*A) + .s^llAf A/f„_i),(l + ^^T''' E ^(^)' ^ « 

fc=n-l /c=0 
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Since the left-hand side is a positive operator, we can write it as B*B. Then, 
1 

Y[{id+sC{ks)){A*A) = {id+sC{ns)){B* B) 

k—n 

n-2 

^ n — 2 

> -s'mU\B*B\\ - s^\\AfMlX\ + ^^T-^ E ^(^)' ■ 

A;=0 

Here, we used p.Sp and the fact that Mt is monotone increasing. This gives the 
following upper bound for ||i?*_B||: 

1 n-2 

\\B*B\\< n ||(id+.£(M)(^*^)ll+5'll^irMf„_i),(l + — ^)"-2E^(^)' 

fc=,i-l fc=0 
1 n-2 



1 " ^ 

fc=n-l fc=0 

1 n-2 

k=n~l k=a 

n-2 

= PIP + .2||A|pA/2^(i + — r-^ 

Tt J ^ ^ 



A;=0 

Therefore we obtain 
1 



W {I + sC{ks)){A* A) 

k—n — 1 

^ n-2 

> -sHi^ArDisr-' - sHii{sHii + i)(i + _^)»-2||^||2 ^ ^(^). 

fc=0 

n-2 

> -s^M^JIAf (1 + -r-'D{sr-' - sHiUl + -T-'\\Af E D{sf 

fc=0 

n — 1 



where to pass to the second inequality we use our assumption on s p.lOp . This 
completes the proof of the bound p.9p . 

To finish the proof of the theorem we use Lemma [T] to approximate the propa- 
gator and put s = in the bound p.9l) . which yields 

^ f hi i"^ 1 ""^ i 

(3.11) n (1 + -^(-))(^*^) > --\\AfMf{i + ^r-i E ^(-)'- 

k—7'L k—0 

Since = (1 + ^^M^f)" ^ 1 as n ^ oo, we get the estimate D{^)'' < 2 for 1 < 

k < n. The right hand side of p.ll|) is bounded from below by ~j^\\A\\'^eMf2n, 
which vanishes in the limit n ^ oo. 
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To show the complete positivity of 74.0 note that any generator CA(t) satisfying 
the assumptions of the theorem can be considered as the generator for a dynamics 
on A ^ B{C^), for any n > 1, which satisfies the same properties, and which 
generates 74.5 ® id acting on ^(2)B(C"). By the arguments given above, these maps 
are positive for all n. Hence, the 74^^ sue completely positive. □ 



4. Lieb-Robinson bound 

Our derivation of the Lieb-Robinson bounds for 7^^ is based on a generalization 
of the strategy |18j for reversible dynamics, and on |25| for irreversible dynamics 
with time- independent generators. This allows us to cover the case of irreversible 
dynamics with time-dependent generators. 

Proof of Theorem [2] Consider the function / : [s, 00) ^ A defined by 

/(i)=/C7t(S), 

where /C G Bx and B G Ay, as in the statement of the theorem. For X C A, let 
X'^ = A \ X and define Cx'= and Cx by 

cx^t) = E ^^(*) 

z,znx=0 
Cx{t) = Cx{t) - Cx^it). 

Clearly, [K., Cx-'it)] = 0. Using this property, we easily derive the following expres- 
sion for the derivative of /: 

- Cx^mjtiB) + lCCx{t)it{B) 

= Cx^{t)f{t)+ICCx{tht(B) , 

Let ■y^g be the cocycle generated by Cx'={t). Then, using the expression for /'(t) 
we find 

m = 7m + I j^:iC£xirhUB)dr ■ 

Since 7^^ is norm-contracting and ||/C||cb is an upper bound for the ||/C|| regarded 
as an operator on A\, for all A, we obtain 

(4.1) ll/WII < ||/(.s)|| + ||/C||eb r \\lx{rhLiB)\\dr. 

J S 

Let us define the quantity 

CB(X,t):^ sup — r-— . 

TGBx II ' llcb 

Note that we use the norm ||7"||cb, because, as mentioned before and in contrast 
to the usual operator norm, it is independent of A. Then, we have the following 
obvious estimate: 

Cb{X,s) < miSviX), 

where Sy{X) = if AT n F = and 5y{X) = 1 otherwise. From the definition of 
the space Bx we get that T{B) = 0, when T G Bx, since B has a support in Y 
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and r n X = 0. 

Therefore ()4.ip implies that 



CB{X,t)<CB{X,s)+ j^\\C-z{s)\\CB{Z,s)ds. 



Iterating this inequahty we find the estimate: 



CB(.X,t)<\\B\\J2 

where: 



n=0 



for n > 1 and ao = 1, (recall that is a constant, that appears in a definition of 
Fp). The following bound immediately follows from this estimate: 



l/Cilc^lS 

Using definition of i^^, we can rewrite this bound as 



11^7^(5)11 < &^||^^||niin(|X|,|r|)e-(''(-'-)-^^(*-^)). 
Op 

So the Lieb-Robinson velocity of the propagation for every t G R is 



□ 



Note that the bound above depends only on the smallest of the supports of the 
two observables. Therefore, in a situation where it makes sense to consider the limit 
of infinite systems, one can get a non-trivial bound when one of the observables 
has finite support but the support of the other is of infinite size (e.g., say half the 
system) . 

We would also like to point out that with the argument given in [22], size of 
the support \X\, can be replaced by a suitable measure of the surface area of the 
support, which gives a better estimate for observables with large supports. 



5. Existence of the thermodynamic limit 

Our proof of existence of the thermodynamic limit mimics the method given in 
the paper [18j . 

Proof of Theorem [3J Denote C„ = Ca„ and 7^^;' ~ 7^"^ Let n > m, then 
Am C A„ since we have the exhausting sequence of subsets in F. We will prove 
that for every observable A G Ax the sequence (7"^(A))„>i is a Cauchy sequence. 
In order to do that for any local observable A 6 Ax we consider the function 

m :=7[':^(A)-7[T'(^)- 
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Calculating the derivative, we obtain 

= £„(t)/(t) + (/:„(t)~/:™(t))7[T^(A). 

The solution to this differential equation is 

fit) = f itXt^nir) - Cra{r)hi"l'^ Adv. 

J s 

Since ^t,r is norm-contracting, from this formula we get the estimate: 
||/WII< f WiC^ir) ~ £,^ir)H":HA)\\dr 



< 



E E rii*zW(7^?(^))iM^- 

-A \ * rir-, ^ J S 



Using the Lieb- Robinson bound and the exponential decay condition p.3p . which 
we assumed holds uniformly in A, we find that 



1I/W1I< 



< 



< 



\\A\ 
\\A\ 



< 



< 



r 



...... E Eii*^wiicbEE^M«^.2/))rf' 

zeAn\A,nZ3z xex yez 

...... j2 EE E ii*^wiicb^M«^:2/))rf^ 

2eA„\A,„ xex yer Z3z,y 



r 



zga„\a,„ xex yer 

3^--'-dr E E^M(rf(^-^)) 

zeA„\A„ xex 

^^C,mkJ'e'^^'--''-dr\X\sup E W^^^))- 

Since F^^ is exponentially decaying when the distance d{x, z) is increasing, we note 
that for n, m oo, the last sum is goes to zero. Thus 

ll(7l"^-7lT^)(^)II^O, asn,m^^. 

Therefore the sequence {lt"s {^)}'n=o Cauchy and hence convergent. Denote the 
limit, and its extension to ^r, as 7]^^. 

To show that 7^ ^ is strongly continuous we notice that for < s < t, r < T, 
and any A G -4}?'^, we have 

for any ri 6 N such that A G Aa„- The strong continuity then follows from the 
strong convergence of 7^"^ to 7]"^, uniformly in s < t G [0,r], and the strong 

continuity of ^[^^ in t. The continuity of the extension of 7]"^ to all of A G At 
follows by the standard density argument. The argument for continuity in the 
second variable, s, is similar. □ 
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